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INTRODUCTION 


The  need  for  a  bubble  radius  measurement  technique  that  is  accurate 

and  does  not  disturb  the  bubble  has  been  brought  about  by  theoretical  and 

1_8 

experimental  studies  of  bubble  oscillations.  For  instance,  theoretical 

work  on  nonlinear  oscillations  of  gas  bubbles  in  liquids  by  Prosperetti* 

predicts  that  individual  bubbles  subjected  to  an  acoustic  field  will  have 

harmonic  resonances  in  their  pulsation  amplitude.  These  resonances  occur 

when  an  oscillating  bubble's  radius  is  equal  to  the  resonance  radius  for 

frequencies  that  are  harmonics  of  the  driving  frequency.  Some 

8 

experimental  work  has  been  presented  by  Crum  and  Prosperetti  to  verify 
these  predictions,  but  the  evidence  is  inconclusive.  A  rise-velocity 
technique  was  used  to  measure  bubble  radii  which,  they  state,  could  be  in 
error  by  as  much  as  10Z. 

Another  example  is  the  prediction  of  the  rectified  diffusion 

threshold  for  various  combinations  of  the  bubble  radius,  oscillation 

frequency,  and  acoustic  pressure  amplitude.*’  Rectified  diffusion  is  the 

process  whereby  a  bubble  will  grow  (under  proper  conditions)  when 

subjected  to  an  acoustic  pressure.  While  there  are  some  indirect 

9 

indications  of  the  validity  of  these  predictions,  direct  measurements  on 
single  bubbles  have  employed  the  rise-velocity  method  which  is  subject  to 
error. 

The  rise-velocity  method  is  the  calculation  of  a  bubble's  radius  from 
a  measurement  of  its  terminal  rise-velocity.  There  are,  however,  certain 


1 


inherent  difficulties  with  this  method 


First,  to  calculate  the  radius 


the  drag  force  acting  on  the  bubble  must  be  known.  The  drag  force  is 
related  to  a  drag  coefficient  and  there  are  dozens  of  different  empirical 
drag  "laws"^  that  are  based  on  rigid  spheres  in  which  compressibility 
and  wall  effects  are  ignored.  As  a  result  large  radius  differences  can 
be  found,  for  the  same  terminal  velocity,  when  using  different  drag  laws. 
Second,  determining  the  terminal  velocity  can  be  difficult.  For  example, 
in  the  method  used  by  Crum,  the  bubble  is  started  from  rest  and  its  time 
to  rise  a  certain  distance  is  recorded.  There  is  a  reaction  time 

involved  in  stopping  the  timer  which  can  be  a  source  of  error  as  great 
as  5%.  Also,  the  bubble's  initial  acceleration  is  assumed  to  be 
negligible.  Third,  bubbles  that  are  large  enough  to  rise  continually 
dissolve.  For  example,  a  20  uto  bubble  is  observed  to  dissolve  in  about 
10  seconds  in  water  chat  is  not  supersaturated  with  gas.  Thus  the 
radius  changes  during  the  measurement.  Accounting  for  this  change  is 
difficult,  since  the  rate  at  which  bubbles  dissolve  is  dependent  on  many 
parameters. **  This  dissolution  is  particularly  a  problem  when  the 
behavior  of  the  bubble  in  a  sound  field  is  being  investigated.  To  let 
the  bubble  rise  "freely"  with  only  buoyancy,  drag,  and  weight  acting  on 
it,  the  sound  field  must  be  turned  off  for  the  period  of  rise.  During 
this  zero  acoustic-pressure  time  the  bubble  will  dissolve  slightly. 

A  different  technique  for  measuring  bubble  radii  is  presented  in  an 

12 

experiment  conducted  by  Marston.  In  it  a  photograph  is  taken  of  the 
light  scattered  from  an  air  bubble  in  water.  This  photo  is  scanned  with 
a  microdensitometer  to  obtain  the  intensity  variation  with  angle.  The 
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photo  pattern  is  then  compared  to  the  pattern  predicted  by  a  physical 
optics  approximation  to  the  Mie  scattering  theory.  This  same  technique 
has  recently  been  refined  and  used  to  correlate  photographic  patterns 
with  predicted  patterns  for  measured  bubble  sizes  from  25  to  1000 
microns. ^ 

It  was  the  difficulties  with  the  rise-velocity  method  coupled  with 
Marston's  work  on  light  scattering  from  air  bubbles  that  led  tq.  the 
present  investigation. 

The  solution  for  scattering  of  electromagnetic  waves  by  a  sphere  is 

generally  referred  to  as  the  "Mie  theory",  named  after  Gustav  Mie, 

14 

although  the  derivation  had  been  worked  out  by  Debye  at  about  the  same 
time.  Mie's  original  paper ^  addresses  the  then  long-standing  problem 
of  accounting  for  the  brilliant  color  exhibited  by  colloidally  dispersed 
metal  particles. 

The  Mie  scattering  theory  is  used  in  various  ways  in  diversified 

1 6-19  20  2 1 

fields  such  as  atmospheric  optics,  astrophysics,  *  — A 


and  physical 


22  23 

chemistry.  Much  of  this  work  is  concerned  with  the  scattering  of 

electromagnetic  radiation  by  distributions  of  particles  of  various  sizes. 

Also,  most  of  the  work  is  concentrated  on  particles  whose  index  of 

refraction  is  greater  than  that  of  their  surroundings.  Only  recently  has 

any  extensive  work  been  carried  out  on  the  case  of  air  bubbles  in 
12  13  24 

water.  ’  *  Studies  on  the  determination  of  the  average  size  of  the 

particles  in  the  distributions  are  based  primarily  on  the  measure  of  the 
angular  intensity  variations.  There  is,  however,  at  least  one  study  in 
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which  Holve  and  Self“  measure  parcicle  size  distributions  using  pulse 


height  analysis  on  near-forward  scattered  light.  Two  excellent  books  have 

been  written  on  the  subject  of  Mie  scattering  and  serve  as  excellent 

26  27 

reviews  of  the  work  previously  conducted.  * 

This  dissertation  deals  with  the  development  of  a  method  to  measure 
the  radius  of  a  single  air  bubble  in  water  through  the  measurement  of 
the  light  intensity  scattered  at  a  single  angle. 

Section  I  contains  the  theory  behind  the  measurements  to  be  made. 
First,  the  equations  which  allow  the  rise-velocity  to  be  used  for  radius 
measurements  are  derived.  Next,  the  theory  for  electromagnetic 
radiation  scattered  by  a  spherical  object  is  briefly  explained.  The 
scattering  functions  are  derived  and  the  far-field  solution  which  is 
needed  to  describe  the  experimental  situation  is  given.  The  calculation 
method  is  described,  and  other  necessary  equations  are  presented. 
Finally,  the  computational  results  are  presented. 

Section  II  presents  the  experimental  portion  of  the  study  and 
consists  of  an  explanation  of  the  apparatus  and  the  procedures  used, 
along  with  the  experimental  results. 

Section  III  discusses  the  results  of  the  experiment.  Some 


suggestions  for  further  studies  and  possible  improvements 
experimental  design  are  also  addressed. 


I.  THEORY 


A.  Rise-velocity  method 


For  a  bubble  rising  in  water  at  terminal  velocity  U,  the  forces  of 

28 

buoyancy,  weight  and  drag  are  balanced.  The  drag  force  is 


i  pU2ua2CD  , 


where  p  is  the  density  of  the  water,  a  is  the  radius  of  the  bubble,  and 
CQ  is  the  drag  coefficient.  The  equation  of  balanced  forces  is 

7ra3og  +  j  pU2ira2CD  *  Tta3pg  ,  (2) 


where  a  is  the  density  of  the  air  in  the  bubble  and  g  is  the 
acceleration  due  to  gravity.  From  this  equation  we  find  that 


a  ,  2  u  c 

a  8  g  CD  * 


since  o/ (p  -  cr)  s  1 . 


The  drag  coefficient  is  expressed  as  a  function  of  the  Reynolds 

number  Re  *  2aU/v ,  where  v  is  the  kinematic  coefficient  of  viscosity 

of  the  fluid  and  is  defined  as  u/p,  where  y  is  the  dynamic  viscosity 

and  P  is  the  density.  The  drag  coefficient  has  been  calculated  by 
28 

Stokes,  for  rigid  spheres,  to  be  24/Re,  which  is  valid  for  low 
Reynolds  numbers.  Many  other  drag  laws  have  been  determined 
empirically.^  The  law  used  for  the  rise-velocity  calculations  in  this 
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study  was  determined  by  Schiller  and  Nauraan  and  was  chosen  because  it 
matched  the  "standard  drag  curve"  fairly  well  for  the  region  of 
interest.  The  form  of  this  law  is 

S  “  fl  (1  +  °-15Re°-687)  •  (4) 

For  comparison  purposes  the  law  determined  by  Langmuir  and  Blodgett^ 

29 

which  was  used  by  Crum  and  Eller  is  given  by 

C-  (l  +  0. 197Re°  * 6  3  +  2.6  *  10“4Re1>38)  .  (5) 

D  Re 

With  each  of  these  drag  laws  an  iterative  method  must  be  used  to 
obtain  the  radius.  The  results  of  such  calculations  for  the  three 
different  drag  laws  given  above  and  a  large  range  of  bubble  sizes  is 
shown  in  Fig.  1.  The  temperature  used  for  these  data  was  25° C.  Note  the 
large  differences  in  radius  for  large  rise  velocities  (about  7  at  0.7 
cm/s  rise  velocity).  This  graph  shows  one  of  the  inherent  problems  in 
using  this  method  to  determine  the  bubble  radius. 


RISE  VELOCITY  (cni/s) 
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B.  Mie  scattering 

Extensive  analysis  of  the  scattering  of  light  by  a  spherical  object 

can  be  found  in  Born  and  Wolfe^  and  Stratton^  as  well  as  in  the 

14  15 

original  works  of  Debye  and  Mie.  The  following  contains  the  essence 
of  these  calculations.  The  potential  approach  will  be  used  rather  than 
dealing  directly  with  the  vector  wave  equation. 

The  system  under  consideration  is  a  non-conducting  sphere 
(medium  1),  radius  a,  immersed  in  a  homogeneous,  isotropic  medium 

(medium  2)  within  which  there  are  no  conduction  currents  nor  free 

charges,  with  a  plane,  linearly  polarized,  monochromatic  wave  impinging 
upon  it.  The  field  equations  under  these  conditions  are 

V  x  E  +  “  0  ,  (6a) 

7  x  ii  -  f  ■  0  ,  (6b) 

V  •  D  -  0  ,  (6c) 

and  V  •  B  -  0  ,  (6d) 


with  6  ■  eE,  and  B  *  uH. 

These  equations  can  be  reduced  to  the  determination  of  a  vector 
potential.  A,  and  a  scalar  potential,  d> ,  from  which  the  fields  are  found 
through  the  use  of  the  equations 


and 


B  -  7  x  A  , 


E 


(7a) 

(7b) 


It  is  also  possible  to  define  an  electromagnetic  field  in  terras  of  a 


single  vector  function.  Assume  that  the  vector  potential  A  is 

proportional  to  the  time  derivative  of  a  vector  IT  1  such  that 


A 


(8) 


Consequently 


B 


ye  V 


(9a) 


and 


E  ■  -  V<$>  -  ye 


_3fii 

3t2  ' 


(9b) 


Substituting  these  expressions  for  E  and  B  into  Eq.  (6b)  results  in 


x  7  x  111  +  7$ 


3  Hu 


(10) 


The  choice  of  the  scalar  potential,  $  ,  is  arbitrary  so  long  as  it 
satisfies 

V2*  -  ye  !^r  -  0  .  (11) 

Therefore,  $  ■  -  7  •  ITi  (12) 

is  chosen.  Integrating  Eq.  (10),  with  respect  to  time,  gives 

3  ^  H 

7  x  7  x  III  “  7(7  •  III)  +  ye  ^ j.2  *  constant.  (13) 

Since  the  particular  value  of  the  constant  does  not  affect  the 


k.  "1  aAeJ*  V*  v’ 


determination  of  the  field,  the  constant  is  chosen  to  be  zero.  Then 


every  solution  to  Eq.  (13)  determines  an  electromagnetic  field  through 
the  equations 


B  -  ye  V  x_^. 


(14a) 


7(7  •  Si)  -  ye 


(14b) 


Using  a  well-known  vector  identity  Eq.  (13)  may  be  rewritten  as 


V2nx  -  ye  -f-y1  -  0 


An  alternate  solution,  ff2 ,  may  also  be  chosen  by  writing 


-U£  V  x  jjp 


(16a) 


7(7  .  n2)  -  ye  f^2 , 


(16b) 


and  by  choosing  the  potentials  to  be 


a  v  * 1  . 

Az  -  ye  — 


3t  ’ 


(17a) 


$2  ■  —  7  •  JI2 


(17b) 


Thus  the  electromagnetic  field  within  our  region  can  be  resolved 
into  two  partial  fields,  one  derived  from  IT L ,  the  other  from  ?i2.  The 
source  of  the  electromagnetic  field  derived  from  the  vector  il x  is  a 
distribution  of  electric  polarization.  This  can  be  seen  by  writing 


D  -  enE  +  where  D  is  the  displacement  current  and  P  the  polarization. 


and  replacing  Eqs.  (6b)  and  (6c)  with 


and 


-  3E  3P 

7  x  H  -  eo  —  -  —  ,  (18) 

7  •  E  -  -  -  7  •  P  .  (19) 


These  two  equations  are  still  identically  satisfied  by  Eq.  (10), 
provided  Hi  is  a  solution  of 


V2jh 


(20) 


In  a  similar  manner,  II z 


can  be  shown  Co  be  a  solution  of 

v2fr 2  -  ye  -|^y2  *  -  m  , 


(21) 


where  M  is  the  magnetization.  This  equation  indicates  that  H2  has  its 
origin  in  magnetic  dipoles. 

Furthermore,  the  vector  potentials  IT  1  and  ff 2  can  be  derived  from  the 
scalar  potentials  n x  and  IT 2  as 

Hi  -  -  7  •  ffx  ,  (22a) 

and  H2  -  -  7  •  ff2  .  (22b) 

Both  Hi  and  n2,  called  the  Hertz-Debye  potentials,  are  solutions  to  the 
scalar  wave  equation 

V2n  -  ye  |^7  -  0  .  (23) 


V  jx.lt. , 


j&L-r. 
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Thus,  the  problem  reduces  to  a  solution  of  Eq.  (23)  in  the  appropriate 
coordinate  system.  For  example,  in  spherical  coordinates  the  fields  can 
be  derived  from^ 


il(rl U)  ^ 


+  k  rlli  , 


f  i  a2(rru)  <2  3 ( rll 2 ) 

9  r  3r90  r  sin0  9$ 

g  =  1  9 2  (rlli)  _  < 2  3  (rfl  2 ) 

<j>  r  sine  9 r 9 <p  r  36 

H  -  32 +  k2rn2  . 


H  <i  9(rni)  1  92(rn2) 

0  r  sin0  3$  r  9r90  ’ 

u  <i  9 Cm i )  1  _  » 

H<j)  “  r  30  r  sin0  3r9$ 


(24a) 


(24b) 


(24c) 


(  24d  ) 


(24e) 


(24f  ) 


where  the  propagation  constant 


with 


k2  «  -  ki<2  “  u2e 


<i  ■  iuie  , 


<2  *  iw 


(25a) 

(25b) 


(25c) 


Since  all  media  are  considered  to  be  nonmagnetic,  u  has  been  dropped. 

For  a  sinusoidal  time  dependence,  e*wt,  the  wave  equation  reduces  to 


72n  +  k2n  -  o  , 


where,  as  before,  k2  ■  oj2e.  In  spherical  coordinates  this  becomes 


The  solution  to  this  equation  is  obtained  in  the  usual  manner  by  the 
method  of  separation  of  variables,  i.e.,  by  assuming 

n  -  R(r)0(e)<t>(4>)  •  (28) 

The  well  known,  separated,  differential  equations  are 


d2rR(r) 


+  [k  - 


n(n  4-  1) 


]  rR(r)  «  0  , 


a hh  lsine  Ir' +  [n<" +  ■>  -  liSei  e<9> 


d2^}~  +  m2$($)  -  0  , 


where  n  is  an  integer  and  m  can  assume  the  values  -n...0...n. 

The  solutions  of  the  radial  equation,  Eq.  (29),  are  the  Ricatti- 
Bessel  functions  defined  as 


n+4(z)  * 

(32) 

n+4(z)  * 

(33) 

where  and  Yn+J^(z)  are  half-integer  order  Bessel  functions  of 

the  first  and  second  kind,  respectively,  and  z  -  kr  (all  function 
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definitions  follow  Abramowitz  and  Stegun  ).  Note  here  that  the  linear 
combination 

tnU)  -  *n(0  -  ixn(z)  -  Hn^!)(z)  ,  (34) 


where  H  ,  '(z)  is  the  half-integer  order  Bessel  function  of  the  third 


kind,  is  also  a  solution  with  the  added  property  of  vanishing  as  z  + 
From  here  on  k^  and  k^  will  be  the  propagation  constants  of  medium  1  and 
2,  respectively,  and  k  is  understood  to  represent  either  k^  or  k^»  The 
solutions  to  Eq.  (30)  are  the  associated  Legendre  polynomials 

0(0)  =  Pn(m)(cos6)  .  (35) 

The  solutions  of  Eq.  (31)  are  the  common  circular  functions  sin(m<f>)  and 
cos(m<t>) . 

The  general  solution  of  the  scalar  wave  equation  in  spherical 
coordinates  may  now  be  obtained  by  a  linear  superposition  of  all  of  the 
particular  solutions,  each  multiplied  by  a  constant  coefficient,  thus 

rn  -  r  l  l  n^m)  ,  (36a) 

n**0  m*-n 

where 

-  [a  ip  (z)  +  b  x  (z)l  (cos0)  [c  cos(m<(>)  +  d  sin(m<|>)]  .  (36b) 

n  n  n  n  n  n  n  n 

Due  to  the  discrete  boundary  at  r  *  a  we  know  that  there  is  not 

only  the  incident  field  (E(i>,  H(i))  and  the  internal  field  (E(w),  H(w)) 

-(a)  -(a) 

but  a  scattered  field  (E  ,  H  )  .  Thus  the  total  electric  field  in 
the  two  regions  is  written  as 

E  ■  E^^  +  E^S^  (outside  the  sphere),  (37a) 

and  E  ■  E^W^  (within  the  sphere),  (37b) 

with  a  similar  expression  for  the  magnetic  field  vector.  Thus  we  must 


also  have  IT  a  '  ' ,  H2 and  Hi'"',  n2'"'  as  well  as  The 
boundary  conditions  require  the  tangential  components  of  the  electric 
field  and  the  normal  component  of  the  magnetic  field  be  continuous  at 
r  ■  a.  Written  in  terms  of  the  potentials,  these  boundary  conditions 
become 


^IrW^W0)!  . 

(38a) 

fj  [r(n2(i)+  n2(s))]  .|j[rn2Wl  . 

(38b) 

im22kor(ir1(l)+  fl2(s))  -  im12k0rH1  <u)  , 

(38c) 

ikBr(I12<l,+  n2<s))  -  ik0rn2(u)  , 

(38d) 

is  the  propagation  constant  in  free  space 

and  mx  and 

m2  are  the  indices  of  refraction  of  medium  1  and  2  respectively. 

To  determine  the  constants  that  satisfy  these  boundary  conditions 
the  potentials  II x  and  H2  must  be  expressed  in  the  form  of  Eq.  (36).  To 
show  that  this  can  be  done,  consider  one  component  of  E  The  radial 
component  of  the  incident  E-field  can  be  written  in  spherical 
coordinates  as 


_  ik2r  cos0  .  „ 

E^  *  e  sm0  cos<t>  . 


Substituting  Eq.  (39)  into  Eq.  (14b)  yields 


ik2r  cos0  .  „  ,  3z(rHi)  .  ,  2_n- 

*  sin0  cos$  =■  — +  k2  rili 


The  first  factor  on  the  left-hand-side  can  be  written,  through  the  use 
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of  Bauer's  formula,  as 


ik2r  cos9  1 


^  l  in(2n  +  1)  ^n(k2r)  P^CcosS) 


ik2r  cos9  .  „  13  ik2r  cos9 

e  ‘  sing  =  -  — - tt  e  4 

ik2r  30 


Using  Eqs.  (41),  (42)  and  the  identities 


P  (cos9)  =  -  (cos9)  , 
30  n  n 


Po(l) (cos9)  =  0  , 


the  left-hand-side  of  Eq.  (40)  becomes 


CO 

eik2r  cos9sin0  cos<£  s  )  £  in  ^(2n  +  1)  <\>  (k2r)  P  1  (cos9)cos$.  (45) 

v*c  2  r  ,  n  n 


This  result  indicates  that  a  solution  of  the  form 


CO 

[i  =  r^T  l  a  \p  (k2r)  P*^  (cos9)cosc 
1C2  »  n  n  n 


is  needed.  Upon  substitution  of  Eqs.  (45)  and  (46)  into  Eq.  (40)  and 
comparing  coefficients  the  relationship 


a  [k2V(k2r)  +  «  in_1(2n  +  1)  ^ 


(k2r) 


is  found.  This  equation  holds  if 


,n-l  2n  +  1  , , 

a  ■  i  -7 — ,  (48) 

n  n(n  +  1) 

since  tp  (kr)  is  itself  a  solution  to  Eq.  (29).  The  calculations  for  the 
n 
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magnetic  potential  are  similar.  Therefore  the  expressions  for  the 
incident  fields  are 

r1ia)  *  ^2  l  1"'1  -nfH"TT  *n(k2r)  P‘‘>(cose,<:OS<’  1 
n=  l 

CO 

and  rn2^  =  -~y  £  iU  *  y-y  ^(kjr)  1  ^  (cos6 )  sin<$  .  (50) 


The  other  potentials  expressed  in  terms  of  a  series  expansion  of  the 
form  of  Eqs.  (49)  and  (50)  with  constant  coefficients  to  be  determined 
are 


rIIi(s)  *  ~  i~r  I  1 


1  r  . n- 1  2n  +  l  .  .  _s  0(0 


(w)  _  1 


n(n  + 

I) 

2n  + 

1 

n  (n  + 

D 

2n  + 

I 

n(n  + 

1)  ' 

2n  + 

1 

n(n  + 

1) 

a  ;  (k.2 r)  P  (cl  s~)cos<J>  .  (51 ) 

n  n  n 


b  ;  (k2r)  P  (cos6)sin<J>  ,  (52) 

n  n  n 


c  .  ( k i r )  P  (cos“)cosi  ,  (53) 

n  n  n 


and  rll 2  ^  =  t~~t  I  i  —7 — --  .-r  d  p  (k2r)  '  (cos9) sinb  .  (54) 

^  n(n  +  1)  rrn  n 


Only  the  function  C  (k2r)  was  used  for  the  scattered  wave  because  it  has 

n 

the  property  of  vanishing  at  infinity.  Also,  only  ^(kir)  was  used  for 

the  internal  wave  because  xn(^ir)  becomes  infinite  at  the  origin. 

Since  the  terms  in  the  series  expansions  are  independent  of  each 

other,  the  boundary  conditions  Eqs.  (38a)-(38d)  must  hold  for  each 

corresponding  term  in  the  series.  This  fact  leads  to  a  set  of  four 

linear  equations  relating  the  coefficients  a  ,b  ,c  and  d  .  These  are 

n  n  n  n 


m[^'(k2a)  -  a  t' (k2a) ]  =  c , 
n  n  n  n  n 


(55a) 


I,  ,.  -:.  -:.  :.  . .  .  .  >  .  .  . 
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C.  Far-field  solution 

Consider  next  the  scattered  field  at  distances  sufficiently  far  from 
the  particle  so  that  k£r  »  (order  of  the  Ricatti-Bessel  functions). 
Under  these  conditions  the  expressions  describing  the  scattered  field 
are  somewhat  simplified^  since  the  Hankel  function  reduces  as  follows: 


?n(k2r)  =  1n+1e-ik!r  , 

(58) 

and 

C'Ckjr)  -  lVik!r  . 

(59) 

A  further  simplification,  in  the  far-field,  results  from  the  observation 
that  the  scattered  wave  becomes  a  transverse  wave  due  to  the  rapid  decay 
of  the  longitudinal  components.  The  transverse  components  of  the  field 
vectors  decay  with  A/r  in  accordance  with  the  inverse  square  dependence 
of  a  spherical  wave  upon  the  radial  distance.  The  radial  components 
fall  off  as  (A/r)2  so  that  they  may  be  neglected  in  the  far-field  zone. 

The  final  result  for  the  far-field  zone  is  then 

H  .  -ikjr 

Ea  -  -  r1  -  —  ■  S2C0S$  ,  (60) 

0  m2  k2r 


H 


and 


E  .  -  — 

9  m2 


6  ie 


-ik2r 


k2r 


Sisin9 


(61) 


In  the  usual  manner  we  identify  the  amplitude  functions 


51  -  l  2,n  x  Iv  (  a  *  (cos0)  +  b  x  (cos0 ) )  , 

L.  n(n  +  1)  ^  n  n  n  n  J 

n*i 

00 

52  =■  V  — 7^-~— rr  (  a  t  (cos@)  +  b  v  (cos0))  , 

i  L ,  n(n  +  1)  v  n  n  n  n 

n=l 


(62) 


and 


(63) 


2< 


where  the  angular  functions  are 


P  ^(cos0) 

•n  (cosB)  =  — : — -  , 

n  sin0 

(64) 

dP  (l)(cos6) 

and 

Tn<cos9)  "  d0 

(65) 

The  energy  flow  in  the  scattered  wave  is 

s  -  7  <EeH*  -  W  '  <66) 

The  intensity  of  the  scattered  radiation  polarized  in  the  9  and  <(> 
directions  is 

1 S2  j  2  cos2$  ,  (67) 

j  S 1 1  2  sin2<f>  .  (68) 

These  components  are  perpendicular  and  parallel,  respectively,  to  the 
scattering  plane.  This  plane  contains  the  incident  direction  and  the 


and 


4ir2r2 


I.  = 


<|>  4ir2r2 


direction  of  the  scattered  wave 


a 


D.  Calculation  methods 


Calculating  the  Mie-scattering  functions  given  in  Eqs.  (67)  and  (68) 
is  not  a  trivial  task.  Many  schemes  have  been  proposed  to  avoid  the 
complexities  of  dealing  with  Bessel  functions  of  large  orders  and 
arguments.  Entire  publications  have  been  devoted  to  presenting  the 


results  of  these  calculations  or  describing  new  algorithms 
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section  will  describe  the  formulas  used  for  the  calculations  of  the 
desired  functions. 

The  desired  quantity  is  the  intensity  of  the  scattered  light, 


I9  -  47*72  |S2|2  cos'* 


h  ■  4^77  lS-l!  sin!5  • 


These  two  components  are  respectively  perpendicular  and  parallel  to  the 
scattering  plane.  For  a  given  wavelength,  in  a  given  scattering  plane, 
and  at  a  constant  distance  away  from  the  sphere,  the  intensity  is 
proportional  to  the  scattering  functions  and  which  are  obtained 


51  *  y  — — —  :t  (  a  f  (cos0)  +  b  x  (cos0)J  , 

L,  n(n  +  1)  n  n  n  n  ' 

n=l 

OO 

52  ■  y  ■>—  X  TT  (  a  T  (cose)  +  b  TT  (cose))  , 

i  L.  n(n  +  1)  v  n  n  n  n  1 

n=  1 


where  the  angular  functions  are 


TT  (COSO)  = 
n 


P  ^  (cos9) 


and 

T  (cos9)  =  — 

n 

n 

d9 

Also  the  coefficients 

a  and  b  are 

found 

from 

n  n 

'P„(a)'P  "(6) 

-  miji 

(3) ^ " (a) 

? (a)t|/'(8)  -  nty  (8)?'(a) 
n  n  n  n 


m^_(a)i^''(g)  -  \p_  (B)\p  '(a) 


n  m?  (a)^'(8)  -  ip  (B)z'(ct) 
n  n  n  a 


where  m  ■  "  ml^m2  *s  t*le  re^at*ve  index  of  refraction,  a  *  k^a  » 
ZTrm^a/Xo,  X0  ■  wavelength  of  the  incident  light  in  vacuum,  B  *  am,  and 
primes  denote  differentiation  with  respect  to  arguments.  The  functions 


ij;  (a)  -  aj  (a)  , 
n  n 

?n(a)  “  ahn(2)(a) 


are  the  Ricatti-Bessel  functions  with  Jn(“)  and  h^  '(a)  the  ordinary 
spherical  Bessel  functions  of  the  first  and  third  kind  respectively. 

By  definition  the  associated  Legendre  polynomial 

dP  (cos6) 

Pn  (cose)  -  -  sine  dcos'§-~  ’  (79> 


Then  by  differentiation 


dP  (l)(cos8) 
n 


=  cose  it  (cose) 
n 


sinz6  tt'(cos0)  =  t  (cos9)  .  (80) 
n  n 


The  recurrence  relationship 
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P  (l)(cos0) 
n 


2n— cosQ  P  (l)(cos0) - IL- r  P  ,(1) 

n  -  1  n-1  n  -  1  n-2 


(cos0)  (31) 


was  used  to  find  IT  by  first  making  a  direct  substitution  to  achieve 
n 


P'(cos0)  -  2n  -  4  cos0  P'  . ( cos 0 ) - -  .  .  - 

n  n-1  n-1  n-1  n-2 


with  Po  ■  0  and  Pi*  1,  and  then  using 


P'  .(cos6)  ,  (82) 


rr  (cos6)  -  -  P'(cos8) 
n  n 


(83) 


Also  by  differentiation 

P;(cose)  -  (2n  -  1)P'  .(cose)  +  P'  _ (cos6)  ,  (84) 

n  n— 1  n-2 

with  Po  *  0  and  Pi  -  0,  which,  through  the  use  of  Eq.  (80),  yields 

Tn(cos0).  These  functions  are  stable  to  upward  recursion. 

The  function  iJj  (a)  was  a  little  more  difficult  to  obtain  since  j_(a) 
n  11 

is  not  stable  to  upward  recursion;  it  grows  rapidly  with  decreasing  n 

(above  the  oscillation  transition).  The  method  adopted  was  a 

40 

combination  of  that  proposed  by  Carbato  and  Uretsky  and 

J.C.P.  Miller.41 

To  avoid  the  difficulties  with  rapidly  increasing  values  during 

downward  recursion  it  was  necessary  to  calculate  the  ratio  r  *  j  ,./j  • 

n  n+l  n 

The  starting  poii;t  for  the  downward  recursion  was  chosen  to  insure 

accuracy  up  to  the  largest  value  of  n  ■  N  needed  for  convergence  of  the 

37 

Mie  series  and  was  found  from 


h- 


\ 

r, 

r, 

r. 


| 

r 

I. 

* 

t 

c 
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N  =  a 


i  / 

+  4.05a  73  +  2  . 


(85) 


Using  this  value  for  N  the  starting  value,  n  -  N',  for  the  downward 
recursion  was  found  from^ 


N'  -  N  -  log,  En[  A  +  %(2.-„y|))l  . 


(86) 


-9 


where  u'  *  2/(2n  +  1),  A  *  0.1,  B  ■  0.35,  and  ^  *  1  x  10  .  Downward 

recursion  was  done  on  the  ratios  using  the  formula 


rn-l  2n  +  1  -  ar  * 

n 


(87) 


with  the  initial  conditions  r„,  -  1  and  r„, -  0.  Recursion  was 

N  N  +l 

continued  until  a  ratio  was  reached,  say  r^,  which  exceeded  unity, 
indicating  that  the  transition  line  had  been  passed.  At  this  point  the 
ratio  was  no  longer  convenient  and  a  change  was  made  to  downward 


recursion  on  a  function  F  using  the  formula 

n 


2n  +  1 


n-1 


F  -  F 


-  r  —  r  .  ,  • 

a  n  n+1 


The  initial  conditions  here  were  F 


X+l 


rx  and  Fx 


(88) 


1.  The  recursion 


was  carried  out  until  n  ■  1.  All  the  F^'s  are  proportional  to  the  real 


j^'s  with  proportionality  constant 


la. 

Fo 


where 


jo 


(89) 


(90) 


sina 

a 


I 


The  real  j^'s  were  Chen  found  by  taking 


j  *  pF 
Jn  n 


up  through  n  ■  X  and  then  using  the  calculated  ratios 


^n+1  rn^n 


for  the  rest  of  the  terms  needed  (through  N'). 

The  above  procedure,  although  somewhat  tedious,  was  the  only 
advantageous  way  the  Bessel  functions  could  be  calculated  for  the  large 
orders  and  arguments  needed. 

The  function  Cn(a)  was  found  with  the  aid  of  the  relationship 


hn  '(<*)  -  Jn(a)  -  iyn(a) 


The  y  ’  s  are  stable  to  upward  recursion  and  can  be  calculated  using 
n 


the  relationship 


Vl(o)  ‘  (2°  ~»~U  yn(c,)  '  yn-l(<l) 


with 


r  ,  ~  -  • 

-l  a 


The  rest  of  the  relationships  used  were 


P  (a)  *  aj  (a)  , 
n  n 

?  (a)  *  ah  ^  (a)  , 

n  n 

^(a)  ■  j'(a)  +  Jn(<*)  » 


*r, 


and 


(100) 


5  (a)  =  ah'  ' (a)  +  h'  ' (a) 
n  n  n 


With 

j '(“)  “  ,  (<*) 

n  n—  i 

n  - 

—  j  (a)  . 

(101) 

QL 

n 

and 

h  (2)(a)  *  hj2)(a) 

n  + 

a 

(102) 

The  Mie-scattering  calculations  were  done  on  an  IBM  4341  computer 
using  the  Fortran  program  which  is  given  in  Appendix  A.  The  program  was 
checked  for  accuracy  by  comparison  to  results  found  in  Refs.  12,  13,  26, 
and  27. 


E.  Computational  results 


For  comparison  purposes  it  was  convenient  to  calculate  an  angular 
efficiency  function  or  angular  "gain"  such  that 


G2 


4  Si 


(103) 


£ 


This  "gain"  is  the  ratio  of  the  scattered  intensity  to  the  intensity 
•  • 

that  would  be  found  in  any  direction  if  the  bubble  scattered  the  entire 

incident  energy  isotropically.  The  parallel  component,  S2.  was  chosen 

for  these  experiments  because  this  component  exhibits  higher  scattering 

intensities,  overall,  and  the  fine  structure  is  less  than  for  the 
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perpendicular  component,  Sj. 

Figures  2a-d  show  the  log(G2)  plotted  as  a  function  of  scattering 
angle,  from  0  to  180  degrees,  for  a  bubble  radius  of  50  pm.  The  graph 
is  presented  in  four  parts  for  added  clarity.  The  parameters  used  in 
these  calculations  are  as  follows:  index  of  refraction  of  the  air  ■  1.0, 
index  of  refraction  of  water  relative  to  air  ■  1.33,  relative  index 

m  *  0.75,  wavelength  Xo  *  632.8  nm.  Figure  2a  shows  the  prominent 
forward  scattering,  sometimes  referred  to  as  the  Mie  effect. 

To  demonstrate  the  sensitivity  of  the  pattern  to  radius,  Fig.  3 
shows  a  plot  of  log(G2)  for  scattering  angles  from  30  to  45  degrees  for 
radii  of  50  pm  and  51  pm.  Figure  4  shows  log(G2)  plotted  for  radii 
between  50  um  and  75  pm  and  a  scattering  angle  of  20  degrees  and  30 
degrees.  This  graph  serves  to  demonstrate  the  need  to  measure 


© 


accurately  the  angle  to  the  detector 


The  choice  of  angles  to  be  used 


Variation  of  gain  with  angle,  radius 


Variation  of  gain  with  angle  radius  =  50ym  and 


BUBBLE  RADIUS  C^rrO 


in  the  experimental  part  of  this  study  was  based  on  a  desire  to  obtain 


as  close  as  possible  a  one-to-one  relationship  between  scattered 
intensity  and  radius.  The  two  angles  decided  upon  have  no  other  special 
significance.  The  two  angles  eventually  chosen  were  55  degrees  and  80 
degrees.  Two  angles  were  selected  so  that  the  55  degree  data  could 
supply  an  initial  calibration  for  the  80  degree  data.  Shown  in  Fig.  5 
is  a  graph  of  the  relative  intensity  as  a  function  of  radius  from  20  ym 
to  100  uni  for  a  scattering  angle  of  55  degrees.  This  angle  was  chosen 
because  of  the  separation  of  the  peaks  and  the  relative  smoothness  of 
the  fine  structure.  Figure  6  shows  a  similar  graph  for  a  scattering 
angle  of  80  degrees.  The  overall  shape  of  this  graph  is  desirable  in 
that  the  intensity  rises  fairly  regularly  (ignoring  fine  structure)  with 
radius. 

In  the  initial  experimental  arrangement  a  large  aperture  (see  Sec. 
II,  Part  A.)  was  used  to  permit  a  rough  trial.  With  this  aperture  in 
place  an  unexpected  shape  was  found  in  the  data  at  55  degrees  (see  Fig. 
14),  in  that  the  low  intensity  points  had  a  larger  amplitude  than  was 
predicted  in  Fig.  5.  In  considering  the  possible  reasons  for  this 
occurrence,  it  was  realized  that  with  a  large  aperture  the  photodiode 
pickup  was  seeing  more  than  just  a  small  angle  around  55  degrees 
(55  ±  .07  degrees).  Subsequent  calculation  of  the  acceptance  angle  for 
this  aperture  (1.6  mm)  gave  2  degrees.  A  numerical  integration  was  then 
conducted  for  scattering  angles  from  54  to  56  degrees.  The  trapezoidal 
method  was  employed  with  a  0.2  degree  step  size.  The  results  of  these 


calculations  are  shown  in  Fig.  7.  The  desired  shape  was  present  with  the 


FIG.  5.  Relative  intensity  as  a  function  of  radius  at  a  scattering  angle  of  55  degrees 


Relative  intensity  as  a  function  of  radius  at  a  scattering 


added  feature  of  reduced  fine  structure.  Similar  calculations  were 
conducted  for  79  to  81  degrees  with  a  step  size  of  0.1  degrees.  The 
results  of  these  calculation  are  shown  in  Fig.  8.  The  almost  total 
elimination  of  fine  structure  was  observed. 

A  functional  dependence  was  needed  to  use  the  light  intensity  at  80 
degrees  as  a  direct  measure  of  the  radius.  A  least  squares  fit  was  done 
on  the  intensity  data  for  radii  from  20  pm  to  100  pm  for  the  integrated 
results  at  80  degrees.  The  functional  dependence  for  this  angle  was 

Radius  -  50 . 1 (Relative  Intensity)*404  ,  (104) 
with  a  goodness  of  fit  of  99.994%. 
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II.  EXPERIMENTAL 


A.  Apparatus 


A  block,  diagram  of  the  experimental  setup  is  shown  in  Fig.  9.  The 
main  components  were  a  7-mW  He-Ne  laser,  a  cylindrical  acoustic 
levitation  cell,  and  a  photodiode.  In  addition  there  were  some  non- 
specialized  items  such  as  the  photodiode  bias  supply  and  amplifier,  and 
the  cell's  driving  oscillator  and  amplifier.  Data  was  taken  and  analysed 
.with  a  Cyborg  Corporation  ISAAC  model  9IA  data  acquisition  system/Apple 
II  computer  combination. 

The  He-Ne  laser  was  a  model  LT-7P  manufactured  by  Aerotech, 
Incorporated,  providing  7-mW  minimum  output  power  at  632.8  nm  in  the 
TEMgg  mode.  There  was  a  500:1  linear  polarization  factor  and  a  1/e2 
beam  width  of  1.0  mm.  The  beam  profile  is  shown  in  Fig.  10.  The  laser 
was  mounted  so  that  the  E-field  was  parallel  to  the  scattering  plane 
(the  scattering  plane  was  chosen  to  be  parallel  to  the  floor)  and  the 
beam  traversed  the  vertical  center  of  the  photodiode  suspension  plate. 

The  cylindrical  cell.  Fig.  11,  consisted  of  two  ceramic  transducers 
separated  by  a  vertical  length  of  glass.  The  cell  was  closed  at  the 
bottom  and  filled  with  filtered,  distilled  water.  The  transducers  were 
driven  in  one  of  their  resonant  modes  at  a  frequency  around  22  kHz. 

The  photodiode  was  an  RCA  device  #  C30957E  with  a  responsivity  of 
0.4  A/W  at  632.8  nm.  The  photodiode  circuit  used  is  shown  in  Fig.  12. 
When  radiant  flux  is  absorbed  in  the  reverse-biased  photodiode,  the 
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FIG.  9.  Block  diagram  of  experimental  setup. 


FIG.  11.  Scale  drawing  (1:1)  of  the  acoustic  levitation  cell. 


electric  field  sweeps  the  created  electron-hole  pair  out  of  the 
depletion  region.  The  resulting  photocurrent  was  used  to  develop  a 
voltage  across  the  external  one  Megaohm  resistor.  This  voltage,  which 
was  proportional  to  the  intensity  of  the  incident  flux,  was  then 
amplified  by  a  high  input  impedance  DC  amplifier,  the  output  of  which 
was  fed  into  the  data  acquisition  system. 

The  photodiode  was  mounted  in  a  custom  built  holder  that  allowed  for 
the  control  of  the  acceptance  angle.  The  holder  was  designed  to  have 
interchangeable  apertures  ranging  down  to  0.1  mm  diameter  for  an 
acceptance  angle  of  0.13  degrees.  The  holder  was  suspended  from  a 
horizontally  mounted  circular  plate.  The  center  of  this  plate  was 
attached  to  a  voltage  divider  circuit.  The  voltage  divider  consisted  of 
a  precision,  linear  potentiometer  which  was  fed  by  a  constant  current 
source.  The  pickoff  voltage  from  this  potentiometer  was  calibrated  to 
provide  angular  information.  The  poi ’ntioraeter-plate  combination  was 
attached  to  a  motor  which  was  used  to  give  angular  sweeps. 


46 


B.  Procedure 

Since  the  scattered  light  intensity  was  critically  dependent  on 
the  scattering  angle  an  accurate  determination  of  the  photodiode  angle 
was  essential.  This  task  was  undertaken  as  the  first  step  in  the 
experiment.  The  zero  scattering  angle  was  defined  as  forward  and  was 
determined  from  a  weighted  average  over  the  beam  profile  of  the  kind 
shown  in  Fig.  10.  The  beam  profile  was  obtained  by  stepping  the 
photodiode,  with  the  motor,  from  one  side  of  the  beam  to  the  other, 
voltage  readings  being  taken  at  each  step.  The  profile  shown  in  Fig.  10 
was  taken  with  the  1.6  mm  aperture  in  place.  Two  independent  methods 
were  then  used  to  determine  the  correlation  between  potentiometer 
voltage  and  scattering  angle.  The  first  method  consisted  of  mounting  a 
spectrometer  at  the  vertical  center  of  the  photodiode  suspension  plate. 
The  telescope  of  the  spectrometer  was  then  used  to  sight  the  angles- 
The  potentiometer  voltage  for  each  angle  was  read  from  a  digital  volt 
meter  and  the  correlation  was  made  by  a  least  squares  fit  from  several 
angles.  The  second  method  consisted  of  mounting  a  diffraction  grating 
on  the  spectrometer,  shining  the  laser  through  the  grating,  and  scanning 
the  photodiode  through  the  pattern.  Matching  the  measured  locations  of 
the  peaks  to  the  angles  given  by  theory  gave  the  relationship  between 
angle  and  voltage.  The  two  methods  agreed  to  within  0.8%  and  provided 
angular  information  accurate  to  within  0.25  degrees.  .4  test  of  the 
relationship  was  made  by  mounting  a  diffraction  slit  in  place  of  the 
grating.  The  results  of  one  of  these  scans  is  presented  in  Fig.  13 


xperimental  points  and  theoretical  curve  of  .04mm  single  slit  diffraction  pattern 


along  with  Che  theoretical  pattern  for  the  given  slit  width. 

The  levitation  cell  was  then  set  in  place  and  aligned  to  be  in  the 
vertical  center  of  the  plate.  A  standing  wave  acoustic  field  was 
established  in  the  cell  such  that  an  antinode  existed  in  the  center. 

i 

This  was  achieved  by  varying  the  water  level  in  the  cell  as  well  as  the 
frequency.  The  radiation  force  on  the  bubble,  in  the  center  of  the 
cell,  was  then  used  to  balance  the  buoyancy  force  and  effectively 
"levitate"  the  bubble,  causing  it  to  remain  relatively  stationary. 
Bubbles  were  initially  induced  in  the  cell  by  increasing  the  pressure 
amplitude  until  gaseous  cavitation  occurred  and  then  reducing  the 
amplitude  to  allow  for  coalescence  into  a  single  bubble.  This  levitation 
technique  also  allowed  for  the  adjustment  of  the  bubble  position  up  and 
down,  through  a  slight  variation  of  the  acoustic-pressure  amplitude. 

The  voltage  corresponding  to  the  light  intensity  was  obtained  as 
follows:  With  the  photodiode  set  at  the  desired  angle  the  position  of 
the  bubble  was  adjusted  until  the  bubble  was  in  the  center  of  the  laser 
beam.  The  ISAAC/Apple  was  then  triggered,  manually,  to  take  the  voltage 
readings  (100  readings/data  point).  The  background  voltage  was  then 
subtracted  from  this  reading. 

The  rise-velocity  radius  was  obtained  by  using  a  sighting  telescope 
that  had  graduated  graticule  marks.  The  actual  distance  that  the  bubble 
traveled  was  known  in  relationship  to  these  marks.  With  the  bubble 
"levitated"  in  the  center  of  the  f ield-of-view  of  the  telescope  the 
sound  field  was  turned  off,  simultaneously,  a  timer  was  turned  on.  The 
bubble  was  allowed  to  rise  a  few  marks,  the  timer  was  stopped  and  the 
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sound  field  turned  back  on.  The  time  and  distance  of  the  rise  gave  the 

29 

rise-velocity,  assuming  negligible  time  to  reach  terminal  velocity. 
The  radius  was  then  determined  through  an  iterative  calculation. 

Little  time  was  allowed  to  elapse  between  the  light  intensity 
reading  and  the  rise-velocity  determination.  This  was  done  in  an 
attempt  to  assure  that  the  radius  did  not  change  significantly  between 
the  two  methods,  since  bubbles  between  20  pm  and  80  pm  change  size  at  a 
rate  on  the  order  of  0.3  pm/s  in  low  amplitude  sound  fields.** 

For  consistency  and  accuracy  certain  procedures  were  always 
followed.  First,  the  photodiode  readings  were  triggered  only  when  the 
bubble  was  in  a  specified  location,  as  viewed  through  the  sighting 
telescope.  Second,  the  rise  time  for  bubbles  was  kept  between  one  and 
two  seconds  whenever  possible.  Times  less  than  one  second  showed 
reaction  time  errors  which  were  noticeable  in  the  large-bubble  data 
scatter  (see  Fig.  16).  Times  greater  than  two  seconds  allowed  the 
bubble  to  dissolve  too  much.  This  effect  was  more  prevalent  in  smaller 
bubbles.  A  20  pm  bubble  was  observed  to  dissolve  in  approximately  10 
seconds.  Third,  attempts  were  made  to  avoid  parallax  errors  in  noting 
rise  distances. 

The  water  was  distilled  and  filtered  through  a  7  pm  filter  to 
eliminate  as  many  background  scatterers  as  possible.  The  water  used  was 
also  slightly  degassed  to  reduce  growth  rates.  The  glass  opposite  the 
photodiode  and  at  zero  degrees  was  covered  (inside  the  cylinder)  with 
black  tape  to  reduce  reflections.  Experiments  were  conducted  with  the 
room  lights  off  and  each  run  was  limited  to  100  data  points  to  reduce 
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C.  Results 

The  photodiode  was  positioned  at  55  degrees  and  rise-velocity  and 
light  scattering  intensity  data  were  taken  (1.6  mm  aperture  used).  The 
raw  values  for  light  intensity  were  divided  by  the  value  which  was 
achieved  for  a  rise-velocity  radius  of  50  um.  The  cumulative  results 
are  given  in  Fig.  14.  These  data  match  that  of  the  theoretical 
predictions  of  Fig.  7  except  for  a  shift  in  radius.  The  spread  of  the 
data  points  was  an  indication  of  the  error  and  error  bars  of  one 
standard  deviation  are  also  presented.  When  the  theoretical  results 
and  the  data  were  compared  it  was  found  that  the  rise-velocity  method 
gave  values  for  the  radius  that  were  10%  too  small.  Figure  15  shows 
this.  In  this  figure  the  value  of  the  rise-velocity  radius  at  51.7  u® 
was  forced  to  match  the  value  of  the  theoretical  curve  at  56.9  pm.  This 
same  10%  shift  in  radii  was  carried  through  for  the  rest  of  the 
experimental  rise-velocity  values.  A  correction  was  also  made  for  the 
shift  in  relative  intensity  due  to  the  initial  error  in  radius.  The 
match  of  the  data  with  the  theory  indicated  that  this  was  indeed  the 
correct  shift.  Since  the  form  of  the  drag  law  is  not  linear  the  actual 
shift  of  all  radii  may  not  be  10%  but  the  overall  fit  indicated  that 
this  was  a  good  first  approximation. 

Similar  measurements  were  made  with  the  photodiode  shifted  to  80 
degrees.  The  same  10%  shift  in  the  rise-velocity  method  radii  was  used; 
the  data  are  plotted  along  with  the  theory  in  Fig.  16.  Again,  the  match 
indicates  that  this  shift  was  real. 


c. 


SCATTERING  ANGLE  -  55  DEGREES 


FIG.  14.  Light  intensity  from  bubbles  as  a  function  of  radius  obtained  by  rise-velocity  method, 
scattering  angle  =  55  degrees. 


III.  DISCUSSION 


The  results  of  this  study  showed  that  accurate  radial  measurements 
of  single  air  bubbles  could  be  achieved  with  a  Mie  scattering  technique. 
Using  the  relative  light  intensity  scattered  at  80  degrees  the 
relationship 

Radius  =  50 . 1 (Relative  Intensity)*404 

could  be  used  to  obtain  the  radius.  This  method  gave  the  radius  of  the 
bubble  to  within  an  estimated  3 %  for  bubbles  of  radii  less  than  80  ym. 

The  estimated  3%  was  a  summation  of  the  possible  sources  of  error  in 
the  light  scattering  measurements.  The  major  source  (as  much  as  2.0 %) 
was  the  ±  0.25  degree  angle  measurement  error.  The  next  largest  source 
was  in  the  voltage  readings.  For  100  readings  the  standard  deviation 
divided  by  the  mean  was  approximately  1.0%,  which  translates  into  a 
less  than  1%  radius  error  for  the  entire  range  of  bubble  sizes.  These 
fluctuations  were  probably  due  to  the  oscillation  of  the  bubble, 
background  scatterers,  changing  bubble  size,  and  noise  in  the 
electronics.  The  numerical  integration  and  the  acceptance  angle 
measurement  added  the  next  most  significant  portion. 

The  other  small  sources  of  error  fell  into  two  catagories.  Those 
due  to  assumptions  made  prior  to  the  calculations  and  those  due  to  the 
experimental  arrangement.  Three  assumptions  were  made  prior  to  the 
computations.  First,  the  index  of  refraction  of  the  water  relative  to 


the  air  was  assumed  to  not  vary  significantly  from  1.33.  Second, 
absorption  of  the  light  by  the  air  in  the  bubble  was  ignored.  Third, 
the  bubble  interface  was  assumed  to  be  pure  air-water,  i.e.  no  surface 
contaminents. 

Some  small  errors  were  inherent  in  the  experimental  setup.  A  small 
component  of  the  perpendicular  field  may  have  been  detected.  This 
component  could  arise  from  inaccurate  horizontal  alignment  of  the  E- 
field.  Inaccurate  leveling  of  the  laser  beam  could  change  the 
scattering  plane  which  would  introduce  a  $  component.  This  problem 
could  also  arise  from  a  deviation  in  the  glass  cylinder.  A  source  of 
error  that  was  probably  most  prevalent  in  the  large  bubbles  (>  80  pm)  is 
the  fact  that  bubbles  oscillating  with  a  great  enough  amplitude  have  an 
average  radius  that  is  larger  than  the  stationary  radius.  1  The  radius 
measured  was  assumed  to  be  the  average  radius,  which  was  equal  to  the 
stationary  radius.  This  effect  may  have  been  a  minor  cause  of  the  data 
scatter  in  the  large  bubble  region. 

Another  source  of  error  was  any  nonlinearity  in  the  photodiode's 
response  curve.  A  rough  check  of  the  photodiode  linearity  was 
conducted.  The  results  did  not  indicate  that  the  response  was  nonlinear 
over  the  range  of  input  intensities  used  in  the  actual  measurements. 
There  was,  however,  a  nonlinearity  in  the  response  over  the  range  used 
for  the  measurements  of  the  beam  profile  and  the  single  slit  diffraction 
pattern.  This  nonlinearity  could  account  for  some  of  the  deviation  of 
the  experimental  points  from  the  theory  in  the  single  slit  diffraction 
pattern  shown  in  Fig.  13.  The  deviation  from  theory  of  this  pattern 


could  also  come  from  the  glass  cover  slips  used  in  the  slit  holder  or 
from  an  error  in  the  slit  width  measurement. 

g 

The  results  presented  here  suggested  that  work  done  by  Crum,  in 
which  the  rise-velocity  method  was  used,  with  a  rigid  sphere  drag  law, 
is  indeed  incorrect  by  as  much  as  10%.  The  causes  of  the  10%  deviation 
of  the  rise-velocity  radius  from  the  light  intensity  radius  in  this 
experiment  could  be  due  to  several  factors.  The  major  source  of  error 
was  considered  to  be  the  choice  of  a  rigid  sphere  drag  law.  The  results 
of  this  experiment  indicated  that  there  was  a  larger  drag  coefficient 
for  air  bubbles,  of  less  than  100  pm,  than  is  predicted  by  a  rigid 
sphere  drag  law.  The  departure  from  rigid  sphere  drag  laws  is  known  to 
exist  for  large  bubbles, ^  but  data  are  absent  for  the  bubble  sizes  used 
in  this  experiment.  Another  source  of  the  deviation  is  the  short  time 
difference  between  light  intensity  and  rise-velocity  measurements.  Also 
any  change  in  the  temperature  of  the  water  during  a  data  run  was 
ignored. 

Although  this  Mie  scattering  method  is  theoretically  feasible  for 

all  bubble  sizes  there  was  an  exceedingly  large  amount  of  data  scatter 

for  large  bubbles.  This  lack  of  correlation  was  assumed  to  be  a  result 

of  experimental  errors  in  the  rise-velocity  measurement  since  the  Mie 

13 

scattering  theory  has  been  shown  to  be  accurate.  The  experimental 

error  was  greatest  for  large  bubbles  due  to  the  fact  that  the  reaction 
time  errors  were  larger  for  large  bubbles  and  the  fact  that  large 
bubbles  were  harder  to  keep  stationary  in  the  laser  beam. 
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Despite  the  lack  of  correlation  for  large  bubbles  the  Mie  scattering 
method  was  superior  to  the  rise-velocity  method  for  measuring  the  radii 
of  air  bubbles,  even  if  the  correct  drag  law  had  been  known.  The  Mie 
scattering  method  was  precise  (estimated  1%  error)  and  once  calibrated 
required  the  measurement  of  only  a  single  parameter.  It  also  did  not 
disturb  the  bubble,  if  light  pressure  was  ignored.  It  was  also  fast; 
data  points  (100  readings/point)  could  be  taken  in  about  0.5  second. 
This  speed  reduces  the  change  in  the  radius  during  a  measurement. 

The  principal  disadvantage  of  this  technique  lies  in  the  need  for 
specialized  equipment  and  precision  optical  alignment.  The  only  other 
drawback  with  the  method  was  the  need  to  use  the  rise-velocity  method 
for  an  initial  calibration.  Before  the  intensity  relative  to  the 
intensity  at  50  um  could  be  calculated  the  intensity  of  50  um  was 
required.  This  was  accomplished  by  making  rise-velocity  radius 
measurements  in  conjunction  with  light  intensity  measurements  at  55 
degrees.  This  55  degree  data  supplied  the  correlation  that  allowed  the 
80  degree  intensity  to  be  used  for  direct  measurements  of  radius. 

Initially,  the  idea  was  to  use  the  ratio  of  the  light  intensity 
readings  from  two  photodiodes,  one  at  55  degrees  the  other  at  80 
degrees.  This  ratio  from  a  single  bubble,  followed  through  a  range  of 
sizes,  would  have  located  the  point  on  the  ratio  curve  accurately.  This 
method  would  also  have  entirely  eliminated  the  need  to  obtain  the  rise- 
velocity.  Unfortunately  a  set  of  matched  photodiodes  were  not  available 
(lack  of  foresight  on  my  part).  With  photodiodes  of  different 
responsivities  it  was  not  possible  to  obtain  ratios  of  one  (for  the  same 


input)  over  the  range  of  inputs. 

Besides  the  use  for  measuring  radius  directly  this  method  has  other 
applications.  With  an  independent  measure  of  the  radius  and  an  accurate 
measure  of  the  rise-velocity,  the  drag  coefficient  can  be  experimentally 
determined.  An  accurate  measure  of  the  radius  of  a  bubble  can  also  be 
used  to  make  a  measurement  of  the  damping  on  an  oscillating  sphere, 
which  has  not  been  accurately  measured.  There  are  two  ways  this  damping 
can  be  measured.  First,  the  radius  at  which  the  harmonic  resonance 
occurs  can  be  found.  This  radius  when  matched  to  the  predicted  radius, 
which  contains  the  damping  and  some  other  measurable  parameters,  will 
give  the  damping.  The  method  may  also  be  refined  to  allow  the  direct 
measurement  of  the  oscillation  amplitude,  which  is  a  function  of  the 
damping.  The  damping  is  of  interest  particularly  for  high  vapor- 
pressure  liquids.^ 

Finally,  there  are  a  few  modifications  which  should  be  done  to  the 
present  setup  to  improve  the  ease,  accuracy,  and  repeatability  of  the 
method.  A  precalibrated  angle  measurement  device  would  reduce  the  error 
by  about  2%.  A  pair  of  matched  photodiodes  and  preamps  are  needed  to 
allow  for  ratio  measurements.  A  wider  laser  beam  would  eliminate  the 
error  caused  by  improperly  centering  the  bubble  in  the  beam.  With  a 
wider  beam  there  also  exists  the  possibility  of  using  a  suspended  glass 
bead  for  calibration  purposes.  Since  the  scattering  pattern  is  sensitive 
to  the  relative  index  of  refraction,  anv  measurements  which  are  to  be 


made  on  other  than  pure  substances  will  require  a  means  of  measuring  the 


index  of  refraction  of  the  mixture.  For  example,  mixtures  of  high  vapor 
pressure  substances  (like  alcohols)  with  water  could  be  used  to  measure 
vapor  effects  on  damping. 
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APPENDIX  A 


Following  is  a  listing  of  the  program  used  to  make  the  Mie 
scattering  calculations.  The  program  when  run  requires  the  input  of 
several  parameters,  M,IRS,IRF,IRI,IDS,IDF,IDI,IOUT. 

M,  is  the  relative  index  of  refraction 

IRS,  is  the  starting  radius  times  10  (ex.  10  microns  is  input  as  100) 
IRF,  is  the  final  radius  times  10 
IRI,  is  the  radius  increments  times  10 

IDS,  is  the  starting  angle  in  degrees  times  100  (ex.  20  degrees  is  input 
as  2000) 

IDF,  is  the  final  angle  times  100 

IDI ,  is  the  degree  increments  times  100 

IOUT,  is  the  output  format  desired,  I0UT*1  gives  deg,G2,  I0UT*2  gives 
radius ,G2 

This  program  was  altered  sligntly  to  give  the  numerical 


integrations 


ill  E 


FORTRAN  A 1  CMS-OLE  MIES 


IMPLICIT  REAL*8  (A-fa,C-Z ) 

IMPLICIT  INTEGER;  I-  N) 

REAL  *8  EJA  '3  000)  ,  RJE  '3  000)  ,E  [3  000)  ,  J  '3000) 

EEAL*8  JO , INTENS , fl 

COMPLEX*  16  Z8,Z9  ,HA0  ,HA  1  ,  R  A2  ,  LIT  A,  LITB,  NWF,TCP,S2,S  1,HPA 
CCP1ELEX*  16  LITA1,LITA2,LITB1,LITP2#SC 
INTEGER  LAM, NEW 
Z8=  »  1.  DO  ,0.  DO) 

Z9  =  ;0-D0,  1.DC) 

PI  =  3.  1  41  59D0 

READ*, M, IRS, IRF ,1  PI ,1 DS , I DF ,IDI , IOUT 
DO  88  IR= IRS , IR  F , 1R I 
PAD=IR*. 1  DO*  1.D-06 
FOUT=IR*. IDO 

A=2.  EO*PI*  1.  33*RAE/f.32  6D-C7 
B=A*M 

BI GN= A+2 .  +  J  4 . 05* ; A**.  333) ) 

N  9=IDI  NT  IBIGN) 

S8=DSIN  ;a)/a 
S9=DC0S (A) /A 
II =1 
X=  A 

GOTO  400 
DO  41  1=1, LAM 
RJ  A  ;i)=PA*J  ;i) 

I  =  L  A  M 

bj  a  ;i+ 1)  =r  ;i)  *r j  a  ;i) 

if  ;rja  ;i+  i)  .lt.  i  .e-77)  goto60 

i=i+i 

GOTO  42 
11=2 
X=E 

GO  10400 
DO  43  1  =  1, LAM 
rje  ;i)  =  pa*j  ;i) 

CONTINUE 
I  =  LAM 

rjb  ;i+  i)  =r  ;i)  *rjb  ;d 

IF  ;RJB(I+ 1)  -  LT.  1.  E-77)  GOTO  70 
1  =  1+  1 
GOTO  44 
CONTINUE 

DO  88  IDE G= IDS, IDF  ,1D1 
S2  =  J0.  DO,  0.  DO) 


© 


S  1  =S  2 

QSCA=0 . ODO 
P0  =  0. odo 
P 1=0. ODO 
F?0=0.  ODO 
P? 1=0. ODO 
Y  A0  =  S8 
YA  1=-S  5 

CEG=IDEG/  100.  DO 
TH=DEG/57.29578D0 

c=  dcos  ;th) 

N  =  1 

48  CONTINUE 

SA=A*RJA  (  N) 

SB=  E*EJB  (N) 

IF  (YA1.LT.-1.0E70) GOT0501 
YA2= (2.DO*N-1.DO)  *YAl/A-YAO 
502  CONTINUE 

HA2=EJA(N)-Z9*YA2 
IF  (N.EC.  1)  GOTO  500 
HA  1=RJA  (N-1)  -Z‘J * Y  A  1 

P2  =  C*P  1*  (  2.  DO*N-  1 .  DO)  /  (N-1  .DO) -PO*N/  (N-1.D0) 
PP 2=  12.DQ*N-1  .DU)  *P  1  +  FFO 
AJP=RJA  (N-1)  —  (N+1.D0)*RJA  (N)/A 
E JP=  RJB  (N—  1 )  —  (N  +  1.DU)  *RJE(N)  /D 
30  CONTINUE 
EIE=-P2 

TA0=-C*P2-  (C  *C-1  .  ODO)  *PP2 
HP  A=HA1-  (N+  1.  DO)  *HA2/A 
SPA=A*AJP+RJA(N) 

SP  B=B*BJP*R JB  (N) 

TCP=  A*  HP  A  +  HA  2 
LI TA 1= (SA  *SPB) -  (M  *S  B*S  FA ) 
LITA2=A*HA2*SPB-M*SE*TCP 
YTESTA=REAL  (LITA2) 

YTESTA=DA5S ( YTESTA) 

IF  (YTESTA.LT.  1.  0E-7  0)  GOTO  503 
50  4  CONTINUE 

LITA=LITA  1/ LIT  A  2 
LIT31=H*SA*SPB-SB*SPA 
LITB  2=M*A  *U A2*SPB-SE*TCP 
YT  ESTB=R  EAL  (LIT  B 2) 

YTESTB=DABS (YTESTB) 

IF  (YTESTB. LT.  1.  E-70)  GOTO505 
506  CONTINUE 

LITE=LITE  1/ LITB  2 
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NWP=LI  TA*TA0+LIT3*PIE 

S2  =  NW  P*  (2 . DO*  N*  1.  CO)  /N/  JN  ♦  1 . DO )  +  S2 

Sl=  (2.D0*N*  1.DO)  /N/  (N*  1.  DC)  *  (LITA*?IE+LITB*TAQ)  +S  1 
SO=£G*  (2.DO*N+ 1.DO)  *  (1ITA  +  LITB) 

TEST=LITA*DCONJG  (LIT  A)  *L  1TB* DCO N JG  (LITB) 

IF  (TEST- LT-  1  -E-2 0)  GCTO50 
CSCA=QSCA+ (2.D0*N  +  1.D0) ♦TEST 
P0  =  P1 
P  1  =P  2 
P?0=PP  1 
F? 1 =PP2 
Y  A  0=Y  A  1 
YA  1  =  YA  2 
H=  N  +  1 

IF  (N.GT.N9)  GOTO  5  0 

GOTO  4 3 

CONTINUE 

h;io=2.*a*  (:i-  1.) 

CSCA=QSCA*2-/A/A 

i:;t2ns=s2*dconjg  ;s2) 

GCNE=S  1*DC0NJG  (S  1) 

G01.E=G0NE*4  ./A/A 
GC  NE=DLOG  1 0  (  GON  E) 

GT  KO=D  LOG  10  ( 4  . *1  N  IE  NS/  A/  A) 

SOUT=SO*  DCON JG (SO) 

SGUT  =  DLOG  10  (SOUT*.  2  5) 

IF  (IOUT.2Q.2)  GOTO  6  00 
WRITE  (6,  10  1)  DEG 
WRITE  (6, 101)  GTWO 
WRITE  (1  1,  100)  DEG,  GONE 
WRITE'  12,  100)  DEG,  GTWC 
GOTO  88 

WRITE  (6,  1  00)  ROUT,  GONE 
WRITE  (11  ,  100)  SOOT,GCNE 
WRITE  ( 12, 100)  ROUT, GTWC 
CONTINUE 
ST  CP 

FORMAT  (2X  ,F6.2,F  1  0.5) 

FORMAT  (2X  ,  F  10.  5) 

FORMAT  (2 X ,3 E  1 3.  6 ) 

FORMAT  (2X  ,14,  2E17.6) 

FORMAT  (2X,I4,E17. 10) 

AJ  F=S8-2.  DQ*RJA  (  1)  /A 

BJP  =  DSIN  (B)  /B-2.  D0*R JB  (1)/3 

HA  1  =  S8-Z9*YA1 

P2=  1  • DO 

P?2  =  0. DO 


►TO 


GOT030 
501  YA2=YA1 
GOTO  502 

400  UPB=2.D0*X/J2.D0*BIGN*1. DO) 

ENEW=DLOG  I1.E-9*  (.1  +  .35*DFF*  [2  . -UPR*  U  PR)  /  2.  /  [  2  . -0  PB*U  PR)  )  ) 
ENEW=EIGN-ENEW/. 693 
NE  K=IDINT  JENEW) 

401  NEV=NEH+25 

402  H(NEW)=0. 

1=  NEW 

4  10  8(1-1)  *X/(2.D0*I*1.D0-X*fi  (I)  ) 

IF  'K  (1-1)  .GT.  1.)  GOTO  420 
IF  (I.EQ.  2)  G OTO  401 
1  =  1-  1 
GOTO  410 
4  20  LA  K=I—  1 

J  JLAM  +  1)  =H  ;LAH) 

J  'L  AF1)  =  1  *  DO 
1  =  LA3 

4  30  J  (I-  1)  =  (2.D0*I+  1  .DO)/X*J  CI)-J  CI+ 

1=1-  1 

IF  'I- HQ-  1)  GOTO  440 
GOTO  430 

440  J0=  C3.D0/X)  *JJ1)  -JC2) 

P A  =  D SI N  ;  X)  /X/JO 
IF  (II. EQ.  1)  GOTO  10 
GOT020 

503  LITA2=Z8*1.D-70+Z9*1. E-66 

GOTO504 

5C5  LITB2=Z8*1.D-70  +  Z9*  I.D-fcfc 

GOTO  506 
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